The low spin moment in LaOFeAs is due to a hidden multipole order caused by spin 

orbital ordering. 
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An antiferro-magnetic (AF) low moment solution, 0.4 /ie/Fe, is found in the case of LaOFeAs 
for an intermediately strong Coulomb interaction U of 2.5-3.0 eV. This solution is stabilized over 
a large moment solution due to the gain in exchange energy in the formation of large multipoles of 
the spin magnetization density. The multipoles are of rank four and can be understood as a kind of 
spin-orbital ordering. Parallels can be drawn to the stabilization of the AF order in e.g. CaCu02- 



With the discovery of the iron pnictide layered su- 
perconductors in 2008 PP, a hope was quickly raised 
that these materials would finally lead to an understand- 
ing of the elusive mechanism of the superconductivity 
of the high-Tc cuprates. Indeed there are many com- 
mon features; the fact that the parent compound is anti- 
ferromagnetic (AF), the central role played by a transi- 
tion metal layer, the fact that the AF order quickly dis- 
appears with doping and then is overtaken by a strong 
superconducting state. However, fairly soon some dif- 
ferences were also discovered. While the main electrons 
in the cuprates are correlated and close to an insulat- 
ing state, in the iron pnictides they seems to be at most 
moderately correlated and metallic [2j [3] . This difference 
between the two types of materials is also manifested by 
the fact that density functional theory (DFT) based cal- 
culations of the undoped iron pnictides obtain the correct 
metallic AF order while in the undoped cuprates they 
falsely lead to a non-magnetic metallic state. This lat- 
ter failure is due to the neglect of the strong correlation 
among the Cu 3d states which is believed to play a crucial 
role in forming the superconducting state in the doped 
materials. For instance, if a correlation term is added to 
the DFT Hamiltonian, local density approximation plus 
added Coulomb U interaction formalism (LDA+C7), an 
AF insulating phase is obtained [4]. However, with the 
ever increasing number of DFT studies, it has been clar- 
ified that DFT has problems also for the iron pnictide 
parent compounds, although of different nature [5 j. The 
calculations systematically overestimate the ordered AF 
spin moment, which is 0.35 /iB in LaOFeAs [6]. In fact, 
state-of-the-art DFT calculations in the generalized gra- 
dient approximation (GGA) give spin moments of the 
order 2.0-2.5 /^b [SI E] ? i-e. an overestimation by at least 
a factor five. 

In this Letter we perform LDA+U calculation for the 
AF parent compound LaOFeAs. The obtained results 
show, that for realistic U parameters, a low spin mo- 
ment solution is stabilized due to polarization of higher 
multipole moments of the spin density. These terms can 
be analyzed as a spin orbital ordering among mainly the 
xz and yz d-orbitals at the Fe sites. Finally we make 
a comparison with the LDA+U solution for an undoped 
cuprate, CaCu02, which reveals a striking similarity in 



the role played by magnetic multipoles. 

The electronic structure is calculated within the 
full-potential augmented plane wave plus local orbital 
(APW+lo) method as implemented in the elk code [8]. 
The LDA+Z7 approach is applied following the same 
methodology as described in Ref. [9] with Yukawa screen- 
ing [10] and around mean-field (AMF) double counting 
while the GGA [TQ is used for the DFT part. The AFM 
Brillouin zone is sampled with 10 x 10 x 6 k points uni- 
formly spaced. The calculations are done for the crystal 
parameters of the experimental high temperature tetrag- 
onal structure [6] . The dimension-less parameter govern- 
ing the number of augmented plane waves R\G + k\ max 
is chosen to be 8.0, where R is the Fe muffin tin radius 
and G are the reciprocal lattice vectors. 

There have been several attempt to estimate the mag- 
nitude of the Coulomb interaction U in this compound. 
The results stretch all the way from fairly large values 
of 4 eV leading to strong correlation [T2] , through mod- 
erate values of 3-4 eV [13 and 2.7 eV [14 , down to less 
than 2 eV [15] . As has been discussed [I3j [14] , part of 
the disagreement stems from the different choices of band 
manifolds that are allowed to interact with this Coulomb 
interaction. If a down-folding down to a subset of Fe 
d-states is performed, the effective Coulomb interaction 
has to be decreased too, otherwise the effect of correla- 
tion is overestimated. In the present study we will vary 
U between and 4 eV, where the eV case corresponds 
to a pure GGA calculation, since all Slater parameters 
are screened with the same Yukawa screening parameter 
[9]. In this approach the Hund's rule exchange param- 
eter J varies automatically between and 1 eV, with 
e.g. J=0.82 eV for [7=2.5 eV, a set of values which is 
very close to the values obtained by a constrained DFT 
approach [14 with [7=2.7 and J=0.79 eV. 

The total energy as a function of the spin moment, 
as obtained by constraining the staggered spin moments 
[T6] of the stripe ordered AF state, and as a function 
of U, is displayed in Fig. [T] In agreement with earlier 
studies [5] the GGA curve (U = 0) has a clear deep min- 
imum at m = 2.2 /iB- This minimum moves slightly to 
larger moments by increasing U. Simultaneously, a sec- 
ond solution starts to develop at a smaller moment. At 
U ~ 2 eV this has evolved to a local minimum, which 



FIG. 1: Total energy per magnetic unit cell (4 formula units) 
as a function of the staggered spin moment per Fe atom cal- 
culated with varying < U < 4 eV in steps of 0.5 eV (some 
values are indicated), with solid curves for integer values. No- 
tice that the energy shifts between the curves are arbitrary 
and chosen such as to simplify the comparison. 



becomes the global minimum for [7 > 2.5 eV, a value 
close to the estimated one [14 . At the largest values of 
the Coulomb parameter also an intermediate minimum 
is formed. The stabilization energy of the low moment 
state is large. Already the high moment solution of GGA 
had a significant stabilization energy of 0.17 eV per for- 
mula unit (fu), for [7=2.5 eV the low moment solution 
is lower than the high moment solution by 0.04 eV. It is 
here worth noticing that except for the pure GGA cal- 
culation the m = solution is not stationary. This is 
an indication that m = generally is not a time reversal 
(TR) symmetric state. In fact it can be much lower in 
energy than the TR-symmetric solution, about 0.2 eV/fu 
for C7=2.5 eV. 

In order to analyze this low moment solution, that is 
stabilized at physical values of the Coulomb parameter [7, 
we will adopt the multipole tensor formalism which has 
been described in some detail earlier [9]. The multipole 
tensors can be obtained from the density matrix through 
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where the matrix elements of the corresponding expan- 
sion matrices T a are given by [9j H8] 
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where a = {kp; qt} is a composite index for the double 
tensor indices k and p and the corresponding components 
q and t, the (. . . )-symbol is the Wigner-3j symbol, and 



FIG. 2: Exchange energy per Fe atom decomposed into multi- 
pole w kp contributions (top, where the number indicates kp) 
and multipole tensor components WqQ (bottom, with num- 
bers indicating q) as a function of the spin moment per atom 
for a fixed [7=2.5 eV. In addition the same quantities are 
shown with symbols obtained by a corresponding calculation 
for CaCuCb with [7=7.0 eV and calculated m=0.59 pb- 



Nkp = nik^sp, where n is the usual normalization factor 
[9j[17]. The operator T transforms the spherical tensor, 
which was used in earlier studies [18], to a tesseral form 
[T9] . In this study we prefer to work with tesseral tensor 
moments w kp since they ensure that the matrices T a are 
all hermitian and hence simpler to interpret. The inter- 
pretation of these multipole tensors are that for even k 
they correspond to the multipoles of the charge (p=0) 
or spin magnetization (p=l), while for odd k they are 
multipoles of the corresponding currents. 

The screened exchange energy takes a very simple and 
appealing form when expressed in terms of these tensor 
moments [9j [18] , 
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Firstly, one note the resemblance with the Stoner ex- 
change, — 7m 2 /4, and indeed one can identify I = Aq 
since m — i^q. Secondly, all other multipole moments 
contribute to the exchange energy in the exact same way; 
with an energy parameter times the multipole moment 
squared. 

The contribution from the exchange energy to the to- 
tal energy curve for [7=2.5 eV of Fig. [l] is decomposed 
in the contribution from the multipoles w kp according 
to Eq. ([3| and displayed in Fig. [2] Besides the spin po- 
larization energy, which is of course quadratic with the 
moment, there is a large exchange contribution from the 
magnetic multipole w 41 . Since it has the largest mag- 
nitude for small moments where it dominates, it is the 
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one that stabilizes the small and intermediate moment 
solutions for large enough U in Fig. [I] The most sig- 
nificant multipole tensor components w a as a function 
of the constrained moment are also displayed in Fig. [I] 
There exist three independent components of w 41 : q = 0, 
2 and 4. They are the symmetry allowed hexadecapoles 
(rank four) of the spin magnetization density which are 
rotationally invariant, two- fold invariant, and four- fold 
invariant, respectively, around a tetragonal axis through 
a Fe site. While q=0 and 4 are both allowed also for a 
local tetragonal symmetry, the q=2 is permissible due to 
lower symmetry at the individual Fe sites caused by the 
striped AF order. These large multipoles result in a very 
anisotropic magnetization density as seen in Fig. [3] for 
the case of £7=3.0 eV, where the magnetization density 
has both large positive and negative values but integrates 
to a small value of 0.3 /ie- 




FIG. 3: Isosurface plots of the magnetization density around 
the Fe sites for the striped AF order and £7=3.0 eV are dis- 
played with positive value indicated with dark/blue and neg- 
ative with light/yellow. The arrows show the directions of the 
small integrated atomic dipole moments. The As atoms situ- 
ated below and above the Fe-plane are displayed with spheres. 

This low moment solution is stabilized through an in- 
tricate competition between gain in screened exchange 
energy and loss in kinetic energy. The gain in exchange 
energy by the Coulomb interaction of Eq.[3]is manifested 
in the orbital dependent exchange potential matrix which 
enters the LDA+Z7 hamiltonian, 
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Again, since Tqq = 1 <g) a z , it is possible to identify the 
Stoner exchange splitting As = Im — AoWqq, that can be 
generalized to the higher multipole splitting A a = AkW a . 
The corresponding T a matrix describes which kind of 
states will split due the multipole a and generally in- 
volves the orbital degrees of freedom. The T matrix for 



the most significant component of the magnetic multi- 
pole, in the orbital basis of xz, yz, xy, x 2 — y 2 and z 2 , is 
given by 
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where a z is the Pauli spin matrix. Hence, the existence 
of moments manifests an ordering of spin-orbit als, 
by e.g. a spin dependent splitting of the xz and yz Fe-d 
orbitals. Similar spin-orbital orderings has been recently 
suggested to play a role in LaOFeAs [20] . 
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FIG. 4: The band structure e(k) along the path X-Y-T in 
the Brillouin zone of the striped AF order (shown in the in- 
set together with the larger non-magnetic BZ with the corre- 
sponding wave vector of the magnetic order M') for varying 
Coulomb parameter U and fixed local Fe moment of 0.4 /is- 
The zero of energy is at the Fermi energy sf • The width of the 
bands are proportional to the expectation value of T^J for the 
band states, with the dark/black (light /red) color specifying 
a negative (positive) value. 

In Fig. [4] the band structures are displayed for three 
different values of U with the local staggered magnetic 
moment constrained to a value representative to low mo- 
ment solution, 0.4 /iB- The /c-path is in the basal plane 
of the Brillouin zone (BZ) starts at X, goes via Y, which 
is half the magnetic ordering wave vector M 7 , and ends 
at T in the center of the BZ (see the inset of Fig. [4|. 
The two points X and Y become equivalent in the case 
of TR symmetry. Firstly, we will take a look at the GGA 
(U=0) solutions. One can see the nesting features of 
the bands at the Fermi energy along the Fr-line. These 
bands originates from the hole Fermi surface originating 
from the T point and the electron Fermi surfaces down- 
folded from the M'-point of the larger TR symmetric BZ 
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[2T] . Surprisingly, these bands seem to be inert to the 
nesting effect, since the bands cross also in the AF struc- 
ture with no hybridization gap opening up. This is due 
to the different orbital character of these nested bands. 

Now, we will focus on the stabilization of the low mo- 
ment solution due to the formation of multipoles. For 
finite U one can see that the bands become polarized 
through the action of the orbital dependent potential of 
Eq. Q, as illustrated by the "fatness" and color of the 
plotted bands in Fig. [4] The width is proportional to 
the magnitude of the expectation value of the matrix T^q 
with light and dark colors indicating positive and nega- 
tive values. Here we can see that there are large splittings 
due to the polarization of the w^q tensor components, 
which results in strong rearrangements, especially along 
the yr-line. This produces a large asymmetry along 
the XY-Yme. This splitting together with the one of the 
iu|q tensor components leads to an effective opening of a 
pseudo-gap at the the Fermi energy and a stabilization 
of this low moment solution. For larger U >4 eV, this 
effect is so strong that a true gap opens up and the so- 
lution becomes insulating. Again, as is clear from the 
intermediate [7=1.5 eV plot, the Fermi surface nesting 
has no direct role in the stabilization of this low moment 
solution, which is in agreement with a recent angular re- 
solved photo-emission spectroscopy experiment [22 . 

Finally, we want underline the fact that the multipole 
needed to stabilize the low moment solution, instead of 
the large moment solution as predicted by GGA, also play 
a crucial role in the formation of an insulating AF solu- 
tion in the cuprates, as e.g. CaCu02- The multipoles and 
their energies for CaCu02 calculated with U =7.0 eV are 
shown in Fig. [2] In this case the existence of the multipole 
is easier to understand as it is essentially a pure x 2 —y 2 or- 
bital that polarizes, giving arise to a non- spherical charge 
and magnetization density. However, the magnitude of 
the multi-poles are of the same order as in LaOFeAs, and 
in fact without these multipoles, the non-magnetic solu- 
tion is more stable. The last is accordance with the fact 
that more exchange energy goes into the formation of the 
multipole than that of the spin moment. Hence, in both 
types of compounds it is the neglect of these multipole 
exchange channels in LDA and GGA that lead to the 
wrong ground state, with either too large (LaOFeAs) or 
too small moments (CaCuC^). This favorable compar- 
ison between the magnetism of the undoped LaOFeAs 
and an undoped cuprate will be explored in further de- 
tails in a future study. Then a crucial issue remains wide 
open; how do these spin and spin-orbital ordered AFM 
ground states of the parent compound, with their sig- 
nificant formation energies, vanish already with a small 
doping, which eventually leads to a high Tq supercon- 
ductivity? One can speculate that the multipole order in 
some form remains beyond the doping where the AFM or- 
der is destroyed, and then constitutes the hidden order of 



the so-called pseudogap region, which is well established 
in the cuprates [23] and has been observed recently for 
the pnictides [24] . 
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